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Abstract 

We generalize the Parikh-Wilczek scheme to the tunneling of a massive charged parti- 
cle from a general spherical charged black hole. We obtain that the tunneling probability 
depends on the energy, the mass and the charge of the particle. In particular, the modified 
Hawking temperature is related to the charge. Only at the leading order approximation can 
the standard Hawking temperature be reproduced. We take the Reissner-Nordstrom black 
hole as an example to clarify our points of view, and find that the accumulation of Hawking 
radiation makes it approach an extreme black hole. 
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1 Introduction 



Whether a black hole can radiate puzzled physicists for a long time after the concept of black 
holes was put forward. In 1974, Hawking proved [1] that black holes are not really black, 
and that they radiate energy continuously and have an emission spectrum of a black body. 
Nowadays, the radiation of black holes is called "Hawking radiation" . Although Hawking 
proposed that the radiation of black holes can be viewed as tunneling, he calculated the 
emission spectrum in light of quantum field theory in curved spacetime without following 
the tunneling picture. 

In 1999, Parikh and Wilczek presented [2] a concise and direct derivation of the Hawking 
radiation as a tunneling process. They suggested that the barrier is created by the tunneling 
particle itself, and utilized the Wentzel-Kramers-Brillouin (WKB) method to compute the 
tunneling rate semiclassically. That is, the WKB expansion was made to the first order, and 
the Hawking temperature was then obtained in terms of equating the tunneling rate to the 
Boltzmann factor. The key point of their scheme is to find out a coordinate system in which 
the metric is well-behaved for calculations at the event horizon. Following the scheme, 
they derived the emission spectrum of massless and uncharged particles from spherically 
symmetric black holes, such as the Schwarzschild black hole and Reissner-Nordstrom black 
hole. Their results show that although the radiation is not precisely thermal, the leading 
order is consistent with the standard Hawking radiation, i.e. the black body spectrum. 
Therefore, Parikh and Wilczek provided a reasonable scheme for the study of the Hawking 
radiation as a tunneling process. However, they just considered the massless and uncharged 
particle's tunneling. Later, an attempt was made [3] for the tunneling of a massive and 
charged particle, in which the main idea is to introduce the concepts of phase velocity and 
group velocity and to equate the radial velocity of the tunneling particle to the phase velocity. 

In addition, based on the Hamilton- Jacobi equation Srinivasan et al. suggested [4] an 
alternative approach to study the Hawking radiation as a tunneling process. As this approach 
is not involved in the trajectory problem, it can be applied to the tunneling particles with 
mass or charge in a natural way. The approach was later developed to extremal and rotating 
black holes [5] , to Taub-NUT black holes [6] and to quantum tunneling beyond semiclassical 
approximation [7]. 

Our motivation is straightforward in this note, i.e. to generalize the Parikh- Wilczek 
scheme [2] to the tunneling process of massive and charged particles from general spherical 
charged black holes. Different from the ways given by ref. [4] and ref. [3], we neither cir- 
cumvent the trajectory problem nor adopt the concepts of phase velocity and group velocity. 
What we want to do is to follow exactly the idea of Parikh and Wilczek and to derive strictly 
the equation of radial motion for the massive charged tunneling particle. We note that in 
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this situation the trajectory of the particle is time-like due to its mass, and is not geodesic 
anymore due to the electromagnetic interaction between the particle and black hole. In par- 
ticular, our result, different from that of ref. [3], can on the one hand recover that of Parikh 
and Wilczek in the limit of zero charge, and on the other hand ensure the dispersion relation 
of the tunneling particle, which coincides with the basic requirement that the observer at 
infinity can detect the particle only on shell. 

This note is organized as follows. In section 2, we derive the equation of radial motion 
for the tunneling particle with mass and charge. In section 3, we calculate the imaginary 
part of the action of the tunneling particle and then give a general expression of the modified 
Hawking temperature. In section 4 we take the Reissner-Nordstrom black hole as an example 
for the application of our generalization of the Parikh- Wilczek scheme. We obtain the 
corresponding tunneling rate and the Hawking temperature, which can clarify explicitly our 
points of view. Finally, we make a conclusion in section 5. 



2 Radial motion of massive charged particle 

A general static spherical charged black hole with mass M and charge Q is depicted by the 
following metric and electromagnetic potential, 

ds 2 = -f(r)dt 2 + g- 1 (r)dr 2 + r 2 dQ 2 , (1) 
A, = h(r)%, (2) 

where f(r), g(r) and h(r) are functions of the parameters: M, Q and r. As was done in 
ref. [2], it is necessary to work in the Painleve coordinate in order to calculate the Hawk- 
ing temperature via particle's tunneling. The Painleve coordinate can be acquired by the 
transformation [8], 

1 1 - n(r) 

dtp = dt + J y ) \ dr, (3) 

V f\ r )w) 

where t p denotes the Painleve time. The resulting Painleve line element and electromagnetic 
potential take the forms, 



ds 2 = -f(r)dt 2 p + dr 2 + 2f(r)J^P-dt p dr + r 2 dn 2 , (4) 

■,/ f(r)g(r) 



The motion of a particle with mass m and charge q in such a spacetime can be described 
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by the following Lagrangian 1 , 



m dx^ dx v dx^ 



(6) 



where r is chosen to be the proper time. Substituting the solution eqs. (4) and (5) into the 
above formula, the Lagrangian can be rewritten as 



m 

£ = 2 



-f{r)t p + r 2 + 2f(r)J 1 , 9 ^ Ur 



f(r)g(r) P 



, 1 — g(r) . 
±qh(r)t p -qh(r) x l j^-^ r, 



(7) 



where the dot stands for the derivative with respect to r, i.e. x M = ^f-. By using the La- 



grangian (eq. (7)) and its corresponding Euler-Lagrange equation, we can derive the following 
equation of motion with respect to coordinate t p , 



~Pt P = = mf(r)t p - mf{r" 



f(r)g(r) 



r — qh(r) = u = const., 



(8) 



where the minus sign before pt p is present due to the positivity of the energy oo of the 
tunneling particle. For a time-like trajectory, we also have 



9lAV^& 1) 



(9) 



which can be explicitly written as 



-f(r)i p 2 + r 2 + 2f(r), 
With eqs. (8) and (10), we can obtain 



1 — g{r) . . ^ 
f{r)g(r) p 



(10) 



1 + qh(r)) ± v 7 ! " W)\l (w + qKr)f - m 2 /(r) 



t, 



mf(r) 



(11) 
(12) 



where the upper (lower) sign in eqs. (11) and (12) corresponds to outgoing (ingoing) trajec- 
tories. As a consequence, we can determine the radial motion with respect to the Painleve 



1 The Lagrangian can be obtained from the covariant equation of motion, 

d 2 x» dx v dxP q dx v 

~d^~ + up ^h~dV~m "~d7' 
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time t p by using eqs. (11) and (12), 



dr ■ -i 

dL ~ p 



\/W)( u + <?M r )) 2 - m 2 g(r) 

= ±/(r) V /( \ t= = ■ (13) 

(u + qh(r)) ± y/l-g(r)y/(u + qh{r)f - m?f(r) 

3 Hawking radiation via particle's tunneling 

To derive the Hawking temperature in the tunneling picture, one has to calculate the imag- 
inary part of the action of the tunneling particle. In ref. [2], this imaginary part for the 
tunneling particle which crosses the horizon outwards from r in to r out is defined as 

pr out rr out rp r pr ut ^jj 

ImS = Im / p r dr — Im / / dp r dr = Im / / -j—dr 

Jr in Jr in JO Jr in Jm df~ 

= -Im / / f dr, (14) 

Jr ln Jm dtp 

where the Hamilton equation jjp = has been utilized, and the minus sign appears due to 
the relation [2]: H = M — u. Note that the tunneling particle is massive in our case, thus the 
lower limit of energy integral is no longer zero but the particle's mass. This consideration 
is indispensable 2 for a massive tunneling particle and we shall see that it will lead to a 
nontrivial effect, i.e. a shift with respect to the energy on the emission spectrum (see the 
last paragraph in section 4 for the details). Moreover, we emphasize that the insertion of the 
lower limit of mass on the integral in eq. (14) does not affect the contour integration 3 . The 
reason is obvious because the contour integration is only associated with the integration of 
dr while the lower limit of mass appears in the integration of du. The concrete processing 
is as follows. We use the Parikh-Wilczek's method [2] which is different from the Hamilton- 
Jacobi method adopted in ref. [9], and in light of the former we exchange the order of the 
integration in eq. (14). When we do the integration of dr, the lower energy integral limit 



2 If the tunneling particle arrives at infinity and is detected by an observer, the trajectory eq. (11) must be 
valid there. For an approximately flat metric, we have /(oo) ~ .9(00) <~ 1, and the electromagnetic potential 
^(00) <~ which leads to h(oo) ~ due to the static electromagnetic source, such as the Reissner- Nordstrom 
black hole we shall discuss in the next section. As a result, we deduce u) > m from eq. (11), which coincides 
with the dispersion relation of the tunneling particle at infinity. 

3 The authors would like to thank the anonymous referee for suggesting a discussion on how the lower 
limit of mass affects the contour integration here, and suggesting an expression for eq. (32) in power of q and 
an explanation about the absence of effects on the deformed Hawking temperature caused by the particle's 
mass in the next section. 
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existing in the integration of du certainly has no effects. Next, when we do the integration 
of du, the lower energy limit then contributes a constant shift as appeared in eq. (19) up 
to the first order of u (see eq. (18)). As a consequence, the integration in eq. (14) is not 
sensitive to the limit m — > 0, i.e. one can take this limit before or after integration. 

The horizon is given by /(r^) = g{rh) = and is a function of M and Q: — 
rh{M, Q). Near the horizon one can expand /(r) and g(r) with respect to the radial coordi- 
nate r and to the horizon r h as the center, 

f(r) = f(r h )(r-r h ) + 0((r-r h )), 

g(r) = g'(r h )(r-r h ) + 0((r-r h )). (15) 
Substituting the above equation into eq. (13), Jr can be approximately expressed as 



dr V f'(r h )g'(r h )(uj + qh(r)) - m 2 f (r h )g'(r h )(r - r h ) 

df~ (r-r h ). (16) 

p (u + qh(r)) + a/1 - g'(r h )(r - r h ) y (u + qh(r)) - m 2 f'(r h ) (r - r h ) 

Considering the self-gravitation of the tunneling particle, eq. (16) should be modified by 
the replacement of M — > M — oj and Q — >■ Q — q. Substituting eq. (16) into eq. (14) and 
switching the order of integration for doing the r integral first, we see that there is a pole at 
r = r h . Completing the r integral by deforming the contour around the pole and using the 
residue method, we have 



Im,S = 7r / -.doj, (17) 

where = r\{M — uj,Q — q) and correspondingly f'(rh) = f ' (r^(M — u,Q — q)) and 
g'i^h) = g' (j'h(M — Cj,Q — q)). In the following, by taking the tunneling probability as 
T ~ e -2imS anc j e q Ua ting it to the Boltzmann factor e~^, we can obtain the Hawking 
temperature from eq. (17). 

Now making the expansions, 

dr h 



f(r h (M-oj,Q-q)) = f'(r h (M,Q-q))-f"(r h (M,Q-q)) 
g'(r h (M-u,Q-q)) = g' (r h (M,Q - q)) - g" (r h (M,Q - q)) 



dM 
dr h 



dM 



ui=0 

Q + 0(Q),(18) 



U!=0 

and substituting them into eq. (17), we see that only the leading order terms, f'(r h (M, Q—q)) 
and g'(rh(M, Q — q)), have contributions to the linear term of u in ImS. As a consequence, 
the Hawking temperature is determined by f'(r h (M,Q — q)) and g'(r h (M,Q — q)) com- 
pletely. Omitting higher order terms in the expansions, we can finish the integral in eq. (17) 
straightforwardly, 

lmS = ^= 2% (uj-m). (19) 

y/f'(r h (M,Q-q))g'(r h (M,Q-q)y 
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Therefore, the tunneling probability takes the form, 



(ui—m) 



T ~ g-SImS 1 _ g y/f'(r h (M,Q- q ))g>(r h (M,Q-q)) _ ^0) 

By equating e~ 2Im5 to the Boltzmann factor e~5% we work out the modified Hawking tem- 
perature, 

~ An • (21) 

Note that eq. (21) in the zeroth order of q recovers the standard Hawking temperature for a 
general static spherical charged black hole. In general, the modified Hawking temperature 
depends on the tunneling particle's charge. 



4 An example: the Reissner-Nordstrom black hole 

For the Reissner-Nordstrom black hole, we have 

fir) = 9{r)=l -m + 9l = (r-r + )( r r.) t 

h(r) = (23) 
r 

where the event horizons r± are 

r± = M± ^M 2 - Q 2 . (24) 

As the tunneling particle crosses the outer horizon, we take r\ = r+. Therefore, the imaginary 
part of the action (eq. (17)) takes the form, 

lmS = n / ^dCj. (25) 

Jm r + -r_ 

Considering the self-gravitation of the tunneling particle with energy Cj and charge q, we 
obtain the event horizons r±, and the initial and final positions of the tunneling particle, r in 
and rout, as follows: 



r± = M-Q± y/(M - Cj) 2 - (Q - q) 2 , (26) 
r m = M+^M 2 -Q 2 , (27) 
r out = M-u+^(M-u) 2 -(Q-q) 2 . (28) 
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Note that 4 eqs. (27) and (28) have been used when we derive eq. (25). Substituting eq. (26) 
into eq. (25) and finishing the integral we have 



ImS* = n 



(M - mf — (M — u) 2 + (M — m)^(M-m) 2 - (Q - q) 2 



-(M - u)^{M-u) 2 -{Q-q) 2 



(29) 



Eq. (29) implies that the emission spectrum is not precisely thermal. As a result, the 
tunneling rate has the form, 

r ~ (T 2ImS ~ e AS ™, (30) 

where ASbh = Sbh(M — u, Q — q) — S B h(M, Q) is the difference of the entropy of the black 
hole after and before the emission. Expanding ImS* in u, q and m to the second order, we 
obtain, 

2 / _ 

qQ (3M 2 - Q 2 " 



ImS = 



TT (M + y/M 2 - Q 2 



1 + 



^M 2 - Q 2 
(M 2 - Q 2 f 2 - MQ 2 



(u — m) 

(m 2 - q 2 ) [m + ^m 2 - g 2 ) 

— 2 (m 2 — u 2 ) + 0(u 2 , m 2 , q 2 , urn, uq, mq) >(31) 



(m 2 - q 2 ) (m + ^m 2 - g 2 ) 

From the first order term of the above equation with respect to u, we work out the corre- 
sponding modified Hawking temperature, 



T H = 



VM 2 - Q 2 



2ir 



(M + a/M 2 - g 2 j 2 + qQ 



3M 2 -Q 2 
M 2 -Q 2 



^m 2 - g 2 



(3M 2 - g 2 )g 



y _ l v ^ g 

2tt (m + v/m 2 - g 2 ) (M 2 - g 2 ) + v/M 2 - Q 2 J 



q + 0{q) 



(32) 



In the following we analyze the characteristics of our result from both the charge and mass 
viewpoints. 

If we set q = in the above equation, which corresponds to a neutral tunneling particle, 
we obtain the standard Hawking temperature. For the Reissner- Nordstrom black hole we give 
the concrete relation between the modified Hawking temperature and the charge of the tun- 
neling particle. This relation implies an interesting inequality, T H \ q Q <0 > T H \ g=0 > T H \ g Q >0 , 



4 It is obvious that we have r in > > r out from eqs. (26), (27) and (28), which gives rise to a pole at 
r = Th in the integrand of eq. (14). Therefore, following the procedure for deriving eq. (17) from eq. (14) we 
obtain eq. (25) for the Reissncr-Nordstrom black hole. 
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which shows that the Reissner-Nordstrom black hole is more likely to emit particles with 
the charge opposite to that of the black hole because a higher temperature corresponds to 
a larger tunneling rate. The accumulation of radiation increases the charge of the black 
hole and finally makes it equal the mass of the black hole. Consequently, the fate of the 
Reissner-Nordstrom black hole is an extreme black hole rather than a Schwarzschild black 
hole if we only consider the effect caused by the Hawking radiation. This is the way that 
the charge of the tunneling particle affects the Reissner-Nordstrom black hole in accordance 
with the modified Hawking temperature, which is novel and different from that of ref. [3]. 
Incidentally, we have made a reasonable assumption \q/Q\ <C 1 when we derive the tempera- 
ture expression eq. (32) from eq. (29), that is, the charge of each tunneling particle should be 
much smaller than that of the black hole. As a result, the temperature eq. (32) is certainly 
positive. 

In addition, it seems to be natural that the particle mass appears in the tunneling rate 
r ~ e~ 2lmS . Because the spacetime tends to the Minkowski type at infinity, the observer 
there detects particles only on shell. As the conserved energy u can be identified with 
the energy of the tunneling particle at infinity, i.e. the dispersion relation is satisfied, u is 
definitely larger than the mass m. This ensures the consistency of eq. (30), which means 
ImS > 0. With eq. (31) in its first order approximation, we see that the mass produces a 
shift with respect to the energy u on the emission spectrum. This is also different from that 
of ref. [3]. We note that it is not an artifact of the expansion that the mass does not correct 
the deformed Hawking temperature. As we can see from eq. (29) which corresponds to the 
expression before being expanded, the introduction of the mass in ImS* has no influence on 
the coefficient of u. (When m = 0, eq. (29) reduces to the formulation associated with the 
tunneling of massless particles.) Precisely speaking, the Hawking temperature is proportional 
to the inverse of the coefficient of u in ImS* (eq. (29)). As a result, the mass of the particle 
has no effects on the modified Hawking temperature. However, it is not trivial that the 
mass produces no effects on the modified Hawking temperature but the shift with respect to 
the energy cu on the emission spectrum. That is, the non-triviality of the effect of the mass 
corresponds to the physical explanation that the particle can be detected only on-shell at 
infinity. 

5 Conclusion 

In this note, we discuss the Hawking radiation of a general spherical charged black hole 
by following the idea of ref. [2] but generalizing it to a quantum tunneling particle with 
mass and charge. The mass of the tunneling particle gives rise to a time-like trajectory, 
and the trajectory is not geodesic because of the electromagnetic interaction between the 
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charged tunneling particle and black hole. Furthermore, when we change the momentum 
integral into the energy integral, see eq. (14), we have to impose a nonzero lower limit of 
energy integral m, see the clarifications under eq. (14). We give the formula of the modified 
Hawking temperature which depends on the charge of the tunneling particle (see eq. (21)). 

Taking the Reissner-Nordstrom black hole as an example, we compute the corresponding 
tunneling rate and modified Hawking temperature, and find that our result is different from 
that of ref. [3]. In our result both the mass and charge of the tunneling particle have effects 
on the emission spectrum which is not precisely thermal either as was first pointed out in 
ref. [2]. To the first order approximation, we deduce a charge- modified Hawking temperature 
(see eq. (32)) which recovers the standard Hawking temperature in the limit of zero charge. 
The modified temperature indicates that the Reissner-Nordstrom black hole approaches an 
extreme black hole under the accumulation of the Hawking radiation effect. As to the effect 
caused by the mass of the tunneling particles, a shift with respect to the energy is produced 
on the emission spectrum. As we know, the observer at infinity can detect the tunneling 
particles only on shell. Consequently, the emission spectrum is expected to have a lower 
cutoff compared with the massless emission spectrum. 

Actually we provide a way to realize the idea of ref. [2] for massive and charged tunneling 
particles. We are considering to apply our method to other kinds of black holes, and shall 
give the results separately. 
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